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ABSTRACT 


This thesis examines a method to determine the depth of 
a point source in an isospeed ocean environment. Using the 
Fourier Transform on the acoustic pressure field in the 
range domain results in the attainment of the acoustic fres- 
sure spectrum in the wavenumber domain and a characteristic 
nodal spacing unigue to the source-receiver depths. 
Quantitative examination of a magnitude plot of the spectrun 
and use of simple mathematical formulae yiela the scurce 
depth. The debilitative effects of narrowband noise and 
surface roughness on the pressure spectrum are aiso éexan- 
ined. The pressure spectrum 1S recognizabie in noise arter 
the pressure field in the range domain nas been lost in the 
noise field. The effect of surface gravity waves on the 
pressure spectrum is similar to that on the pressure field 
in the range domain: the characteristic nodal spacing is 


Suppressed as the height of the surface waves increases. 
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This thesis is the third in a series of investigations 
into a propoSal put forth by R. Lauer of Naval Underwater 
Systems Center, New London, CT, in 1979. In his memorandun 
[Ref. 1], Lauer describes a "new" way to analyze sound 
propagation in the ocean and two applications of the methcd, 
which he named the "Wavenumber Technique," or f.T. In his 
proposal, Lauer stated that a source of sound in the ccean 
could be pinpointed in both depth and range Dy a Ssingie 
omnidirectional hydrophone provided the source is generating 
a continuous wave tone. 

The W.T. waS initially described by Fr. DiNapoli [Ref. 2] 
aS an intermediary step in his development of a speedier 
computer algorithm used in analyzing sound propagation as a 
function of Lange. DiNapoli and Lauer proposed converting 
the pressure aS a function Of Tange Pit) to tie beso ae 
Spectrum aS a function of the wavenumber P(k). She cobyver:— 
Sion is accomplished by taking a weighted Fourier transfors 
or the pressure mabe ee P8118) Once the pressure Spectrum is 
obtained, analysis of the pressure field can be accomplished 
in a manner analogous to that presently used in signal 
processing. 

Lauer proposed two uses for his "Wavenumber Techniques" 

1. determination of the depth of the acoustic Source, an 
ability which has obvious tactical applications; 

2. use of the W.1. to evaluate the accuracy of existing, 
and future, acoustic models such as P.E. and 
F.A.C.7T., by breaking the generated pressure field 
down into ray path "families," such aS bottom-bounce, 
refracted-surface-reflected, and surface-ducted fami- 


lies, thus giving a quantitative read-out of what 
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proportion of the totai acoustic energy is being 
channelled through the various ray paths. 
Hopemoth uses, a knowledge of the acoustic environment is 
necesSary. 

It 1s the intended purpose of this thesis to investigate 
the validity of the first application of Lauer's wavenunber 
feemnrgue by Susing it im Meonjunction with an acoustic 
phenomenon that is well-known and for which accurate results 
Ganmepe calculated with orecision; this phenomenon 1s the 
Peeyas Mirror for an Acoustic Doublet. Using several 
source/receiver combinations, this investigator intends to 
eoivpaleweactual Somece depth wath that predicted by the 4.T. 
Tf Lauer's tecanigue is cogent, the source depths calculated 
via the two methods should be equal, or very close within a 
statisticaliv acceptabie degree. 


AS mentioned in the beginning of tais thesis, two prior 


anvestigators, 6B. Stamey and Je Blanchar:, looxed into 
applying the W.T. to determining source depth. Stamney's 
investigation (Reg. 3) utilized Aeeeoe Ina DOUG ecuation 


computer model developed by H.K. Brock, the Split Step fast 
Fourier Transform, or  SSFFT, to generate the pressure 
fields. He ygualitatively compared the model to a Normal 
Mode model and a P. E. Finite Difference model, both range 
dependent, and to DiNapoli's Fast Field Program model, a 
Tange independent aodel. Stamey concluded that the W.T. 
Showed promise aS an analysis tool and that further investi- 
gations were warranted. 

J. Blanchard carried Stamey'S researches one step 
further [Ref. 4]. Using two Parabolic Equation computer 
elcortenns, Brock*’s Spitt Step Fast Fourier’ Transform anc 
Jaeger's Implicit Finite Difference, as his pressure field 
generators, he examined the use of the spaciny between nulis 
of the pressure spectrum to determine source depth. His 
results are interesting and support the need for further 


investigation. 
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The author Wa S Snot Satisfied with Stamey's and 
Blanchard's £indangs. Both individuals acknowledged short- 
comings in their respective studies, especially with regard 
to the mutual presence of the "U-shaped phenomenon" encoun- 
tered in the pressure spectrum. One should remember, 
however, that their theses were preliminary studies only and 
were produced within a highly restricted time frame. Also, 
they used existing ccmputer algorithms specificaily designed 
to approximate acoustic propagation in am velocity—-variant 
mediun. However, the environmental settings used by both 
gentlemen were that of the Lloyds Marror for an Acoustic 
Doublet which reguires a constant sound speed throughout the 
water column; thus, the approximations, assumptions and 
"fudge factors" used by these models make their results for 
an isospeed medium highly artificial and somewhat susrect. 
Fortunately, however, there exists for the Lloyds Mirror a 
Simple, geometric soiution specifically designed for an 
isospeed environment [{Ref. 5]. Since the pressure fields 
for varying source and receiver depths and, therefore, the 
corresponding pressure spectra, can be precisely caiculated, 
it was thought this model would be a good check on the oper- 
ational applicability of Lauer's proposed techniyue. 

Secondary considerations of this study were to examine 
what effect, 1f any, the introductions of, firstly, surface 
waves and, secondly, noise would have on the pressure spec- 
trum. The first objective was Simulated by use of 
Mathematical formulae given in reference 5; for the second 
objective, it was thought that complex noise, random in both 
amplitude and phase and Similar to a Rician distribution 
{Ref. Gov) p.) c sae would provide a simple but reasonable 
approximation of an ccean noise field. This investigator 
intends to gradually intensify the surface waves and the 
noise fieid, separately, until the original pressure spec- 


trum is no longer recognizable. In this manner, it shouic 
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be possible to make a qualitative assessment of their 
respective debilitating effects under the carefully 
controlled conditions found in the Lloyds Mirror phenomenon. 

Section IT sets forth in more detail the theoretical 
development of each point described in this introduction. 
Section III summarizes the investigator's results and the 
conciusions drawn frcem those results. Section IV contains a 


listing of the computer algorithms utilized. 


Vai 


TY. THEORY 


This chapter presents the mathematical basis of the 


wavenumber tecanigque. 


Aw. THE LLOYDS ATRROR PHENOMENON 


A diverging mcnorrequency spherical pressure wave 


[Ref. 5: p.112] can be written in conplex form as 


i(kR- Gt) 


j A 
PIR fi = lire te ie (eqn 2.1) 


where P(K) 15 the spacial factor. Referring to Figuee Zea 


R, 18 the fEange from the Ssounce to Ulcer ec cein ie ea ees 


2 
Lad 
oe 


the range from the image to the receiver. ror convenienc 

in ali that follows, we saall Set A to Unie wagr eevee. The 
ng 

directly out irom the source PS and 22om Uheeapparciue 


following ecuations apply to the acoustic waves propagati 


ee (which is actually the wave reflected from the surface): 


2; / (kR,- WEI 
BLR,t) = L(R)et' = Re | 


DORA = PIR) e 1%! = pa” (eqn 2.2) 


where PR o« Bs » and Fi are compiex functions of horizontal 
Tange £,vertical depth Z and time €; “tho minusss cee 
equation for P; is derived from the surface reflection coef- 
ficient, joy = ~1, EOE ‘a Sincere cuetace Inspection of 
Figure 2.1 reveals that 


/ a 2 
Hy (Zee ré Ro= V(Z-+2Z5) +6 y 


> 
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mance, Af = | R- Rial 


Prose the total field can bDeewEl_tten as 


CkR -Wt) i({kKR,- Wt) 
-iWt e! e . 
Ps ‘ As C R, R, (eg ieeZ. 2 ) 


and oy both have the same time factor, 


Since pop I 
~S ~ 
exp (-iwt), we can retain just the spatial factors and equa- 
(roms 2.3 reduce to 


mi) Ci a 


R, R. 
(eqn 2.4) 


tion 2.4 1s a form of the complex pressure as a function 


tr 


Cat 2 he 
of range used in the computer aigoritha shown in Appendix Tf. 
HOES R22 | tes line 


rh 


Tnspection of Figure 2.1 reveals tha 


a very smail, 


Sing 
GE 
Siig le: ~ a 
or 
Ar = 42s 
R (ecn 2.5) 


and the pressure aaplitude can be approximated by 


in(Zezs 
(eqn 226) 
for the pressure amplitude, 


Looking at just the formula 


we can see that as 
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KZ, 
“ee SS oe J N=O,i¥2e,9 eo 6 (equ Zot) 


the pressure amplitude goes to zero, producing the classic 


JP(K) |] vs R curve shown in Figure 2.2 


B. TBE RELATIONSHIP EETWEEN K, VY, AND BL 


in his Gescripticn of the Wavenumber Techniygue [Ref. 1: 
p- 5-6), Lauer utilizes the horizontal and vertical compo- 
nents of tae wavenumber, k. The general relationship amoung 
these three terms is illustrated in Figure 2.3 and can be 


written nathematically as 


k = the wavenumber = 27f/c 
'Y = the horizontal component of k = k cosY) 
BL = the vertical component of k = k sin D 
Or 
ko 7¥*+0° (eqn 2.8) 


C. THE WAVENUMBER TECHNIQUE AND THE LLOYDS MIRROR 


Given a point source in free space, the monozregquency 
pressure field at the receiver, P(R) , can be expressed as a 


spherical wave (time factor suppressed), 
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PCRS see (ecn 2.9) 


Moere, Lromn Figure 2.1, 


or 
(eqn 2. 19) 


in integral BO eee) ed io wwii Gon ewe heh 7 > epee 7 | 





PIR) = 





> Gacy 


@IkR Cc J.(¥6) et L12r- 2s] 
a 3 if (eqn 2.11) 


Weteg 1S taken from the Fourier-Bessel Transtorm Pair 
(rer. 7; pe 126), 


£0) = [Paar Vd 7 


gins [Frid rr rar fegn 2. 12) 
O 


Miser) ssepresents the acoltistic pressure function in the 
range domain, and BLE) is the acoustic pressure spectrum in 
the wavenumber domain. The sign of the exponential function 
wedestntegral 15 based on which of Ze and Z, is greater. 
For the case of receiver depth being greater than source 
depth so that waves from both image and source are travel- 
ling downward at the receiver depth, the total pressure at 


the receiver can be expressed as 


02 iBl2zp-Ze5) _ eiBl2r + Zz. 


) 
piay= | Len yay 
O I f3 


(egn Z. 13) 
Msewot Euler's identity reduces Equation 2.13 to 
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CO Fi 2 ee 
PIR) = -2 [ et sin(B2Z~) oo wiae ee (eqn 2. 14) 
O 


Usiny the relationship between the Hankel functions [Ref. 5: 
Pp. 449], 


He(yr) = d,lyr) + i¥yr) 
and 


Regie 
Alyn) = Llyn - Yer) 
(eam Z. Io) 


one Can re-write the Eessel function in equation 2.1714 as 


(I) 
J,(yr) = = EP OWe) os HY tyr) | (eqn 2.16) 


so that 


me coy (1) oo (2) 
2qty) = [fH Grar + [BOM trIrdr oan 2.47) 


Letting rt = -r and looking at the {second tegn ont iewsicat 


hand side of eyguation 2.17, 


Ss RE = P 5 : | 
rf flr’) He (yr) r’dr = [fC0KS tyr ar (eqn 2. 18) 


Now, assuming that fir) = E(-r) , then 


_ nes (1) 
ZUG) = [ £i A, (yr) r dr (eqn 2. 19) 
Thererore, 
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(/) 
Dip) = — ay H, (yr) 2/B2r : 
w8 ae sin(GZ,) dy cree 


EOE yr > 2), the aSymptotic approximation of the 
Fankel function for large argument can be used: 


(1) (yr) 
lab ( ie WYP -"g 
Va Te (eqn 2.21) 


aSsubing |c i(r) | goes to zero Paster than (Ln a, 


co 


PIR) =- /—*%- e'4 asa eM, Si O4Ay) el M gy (eqn 2.22) 
Ca iG) s 


-€ 


bet the terms inside the brackets in equation 2.22 be 
ye =~ , , 
defined as aly). then 


His {ag ta oe | 
EIR) = — ia tT] gtyelTdy 


and 





if giyle'Vdy = -4] Meevr LCR) (eqn 2.23) 


Gepeisrng equation 2.12, 


Gag = — 4 (x Poa f VPS (EN e e Maa 


(econ 2.24) 


teee tet eduation 2.24 has the form of the Fourter trans- 


form cf the pressure function if we define 


ite = yer (egn 2.25) 
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SO thatr/r P({R) and wy) are Fourier Transform pairs. 


This equation can be easily evaluated [Ref. 8], 


= 7 | 
NOD) = he. amet sin(224)| 


(egn 2.26) 


HOw Zp > Zc ~- And the magnitude of aty) is 





Si) (ioe 


cae 
| 9(7)) 2 (Cai eZ we) 


For the case where 2, > Zy ; it can be shown taat 
site 
i* WE ibe 
7) = — “ETB sin(Gz J. a 


and the magnitude of ary) is 
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a(n|= 


A ciose look at the magnitude cf the pressure spectrum func-~ 





feqcn 2.29) 


tion in equation 2.27 reveals that nodes occur for values of 


BZ = 07. , 2 S02 ee (egn 2.30) 
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If the magnitude of the pressure spectrum is plotted asa 
runction of the vertical wavenumber then the spacing between 
nodes AG is uniform and the depth of the source can be 
derived from the relationship given by 


—> 


Spee 
“s ~ AG (eqn 2. 31) 


In the event the source is deeper than the receiver, looking 
at equation 2.29, we can see that nodes now occur for values 
of 


BZ = ni ) gear 0 Se. 9/4 aera (eqn ye 3 2) 


Plotting equation 2.29 as a function of the vertical waven- 
unber, the spacing between nedes will now reveal the 


receiver depth based on 





Sa Aa (egn 2. 33) 


ieee 10 Equations 2.31 and 2.33 the exchange of source and 
receiver depths. By placing the receiver at a sShailower 
depth than tne source, no new information (namely the scurce 
depth) is to be found. 

In summary, there is away to calculate the complex 
pressure amplitude as a function of range uSing the complex 
pressure spectrum. Ccnversely, if the complex pressure P(r) 
has been measured at the receiver, then the conplex pressure 
spectrum P(k) can be derived by making use of the relation- 
Ship between the Fourier-Bessel Transform Pairs (see equa- 
tion 2.12), and the depth of the source can he found froa 


the magnitude of that spectrum (see equation 2.31). 
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The foregoing is the mathematics involved in deriving 
the theoreticai pressure spectrun. One purpose of this 
thesis was to compare the theoretical value of the pressure 
spectrum as derived with the Fourier-Bessel Transtcrm with 
the pressure spectrum obtained by use of the computerized 
Discrete Fourier Transform (DFT), otherwise known as the 
Fast Fourier Transform, or FFT. At this Deane, it will be 
helpful to review what happens in the FFT. 

Starting with the Fourier-Bessel transform pair shown in 
equation ZaiZ i (Lr) and gly) can be re-written in terms 


of a discrete sum as is done in the DFT: 


ae ey j2imoa 
FFT “\ gly)} = y ginigare: = £(mAr) 
n=0 {eqn 2. 34) 
where N = 1,2,3,..., up to some large poSitive integer, and 


represents the humber of teanstiorm ports ane 


emmn 


Y f= er ey (eqn 2285)) 


based on the relationship between Ac and AY as described 
in sampling theory [ Retse2-p. Zi, 


OF 
oT (eqn "2. 36) 


By convention, the sign of the exponential function in equa- 
tion 2.34 is taken as negative when performing the "forward" 


tranrsionny, 


PPT Lima SS SGuneme (eqn 2.37) 
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and is taken as positive when performing the "inverse" 


EranmsStor ni, 


FFT“ ainay) = f(mAr) 
(CqneZaoe) 


Then, 





nh Gipe — 
gin hy) a cade (G25) 


Dmrerefore 











elle 
a) = —— e '# ay FFT Oa) | 
and | J 
| gindy)| = a SHG Gira) (eqn 2.39) 


D. THE EFFECTS OF SURFACE ROUGHNESS 


Surface roughness [Ref. 5: p.409] reduces the Lioyds 
minor erfect. The rougher the surface, the greater the 
effect for angles increasingiy closer to grazing incidence. 
Iii M represents the surface roughness factor, then equation 


2.3 Can be re-written as 








ikR 2 _ikR 
-M 
en a i 
where 
M = a sine. (eqn 2.406) 
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PD is the grazing angle of incidence, Bois the average 
height of the surface gravity waves, and IN is the 
wavelength of the narrowband continuouS wave acoustic 
Signal. When M is less than 1, the surface 1s said to be 
smooth; when M is greater than 1, the surface is said to be 
BoOwGo . 

Looking at Figure 2.1, we can see that the overall 
effect of increasing surface roughness is to reduce the 
contribution of the surface reflection or "image" to the 
interference pattern at the receiver. In other words, as H 
increases, the energy detected by the receiver is increas- 
ingly only the energy coming directly from the source. Wher 
the effect of surface roughness 1S a maxinum such that the 
contribution to the pressure field of the “image" is 


completely suppressed, then 


co ten 
gy) = | 2 lyr) rdr 
! 


2 iBl2s -2r| 
= (eqn 2.41) 


me; 


and the magnitude of the spectrum is 





] 


B (eqn 2. 42) 


MGs = 


Ee. THE EFFECT OF ADDING NOISE 


In an effort to simuiate a realistic ocean environnent, 
a normally distributed narrowband noise field based cn a 
Rician distribution [Ref. 6: p. 189} was adapted. The 


resulting acoustic pressure field can be expressed as 
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wei ime) SL CMEY. 8g) (egn 2.43) 


where fA is a hormally-distributed random function with a 
mean ci zero, a standard deviation of one and possesses both 
amplitude and phase. The noise operates independently on 
both the amplitude and phase components of the pressure 
fieid. Since n is independent of range and wavenumber, it 
is treated as a constant by the Fourier Transforn. 

A fourth intention of this thesis was to discover the 
degree to which ncise degrades the waverumber spectrun 
compared with the corresponding degradation of the range 
domain. To produce the random noise fields used in this 
research, the routine listed as "GGNML" in the IMSL likrary 
[Ref. 9] was called upon twice to generate independent, 
pseudo-random functions which interact separately and sinul- 
taneously with both the amplitude and phase components o£ 
P{R) z Whate NOt an elaborate Scheme, owing to time anda 
resource constraints, it was felt that this modest sinu- 
lacrum of narrowband noise would give a fairiy accurate, 
first cut "feel" for the effects of noise on the wavenumber 
Spectrum. The reader should be aware, howevel, that no 
direct consideration cf coherency along the range path was 


taken into account by this method. 
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The Geometry of the Lloyds Mirror Effect 


Figure 2.1 
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Prgures2. 2 


Bee, che FRI ALGORITHSA 


The research for this thesis was conducted entirely with 
computer algorithms to model the Lloyds Mirror, the Fast 
Fourier Transforn, the range windows, waves, and the narrow- 
band noise. This was necessary since the Lloyds Mirror is 
an idealized representation of a situation that occurs onlv 
rarely in nature, and at that is linited to isospeed surface 
ducts. Also, time and money constraints were such that the 
use of computer modeis was an absolute necessity. For these 
reaseons and because the FFT is the heart of the Wavenunber 
Eeemuicwe, — finding a feliable, easy-to-use computer algo- 
rithm to perforu the FFT on the complex pressure function 
was deemed very important. 

inset pe the IMSL libcary routines, FFT2C and FFTCC 
EReE. 93 ome Pg were used to generate the pressure 
spectra. Although the author was unable to ascertain just 
Ney Stamey performed the FFT on his data, it is highly prob- 
deo thats ne used an IMSL iibrary routine. APA ey dees 
thesis, Blanchard stated that he used the FrFTCC routine. 
The magnitude of a typical spectrum generated hy the routine 
FEPZC LS ShOwn in Figure 3. 1. 

The same scenario used to generate Figure 3.1 was used 
EOupzoauce the Spectrum shown in Figure 3.2, but the routine 
PElGG Was Utiliazed. AS Can be seen, the resulting graghs of 
the pressure spectrum magnitude are practically identical. 
Indeed, the set-up of the pressure fields for insertion into 
Saem GOutine differ  Cnly im that FFT2C requires a data set 
COMSiSting @£ an integral power of two number of points, 


whereas FFTCC will handle any number of points. Anc, when 


ES: 


the Cooley-~Tukey FFT listed in Appendix I was used to trans- 
form the same pressure field as was utilized to produce 
Figures 3.1 and 3.2, Figure 3.3 was the result. This algo- 
rithm differs somewhat from the IMSL routines by the manner 
in which the pressure field data is treated. Both FFT2C and 
FFTCC are designed to handle a complex array. The 
Cooley-Tukey FFT, on the other hand, is designed to trans- 
form real data; but the pressure fieid is complex. To get 
around this apparent conundrum, the pressure data is defined 
as either real or imaginary by its placement in the array. 
AS can be seen, the resultant graph of the spectrum magni- 
tude is almost indistinguishable from that produced by 
either IMSL routine. in -eonclusren, 1t does not matter 
-which FFT routine is used to generate the pressure spectrum, 
as long as the pressure field in the range domain is 
"fitted" into the array in a manner suitable to the 


particular algorithm used. 


TABLE I 
Critical Values Used in the Research 





K (m7+y N r_(m) (m7) 
3.0 1024 0.524 1.172 x 1072 
280 1024 0.785 7.816 x 107° 
1.0 ODE lise 3.906 x 107? 


It 1S important to note four Chings Vaboutstno 9 Fpaae 
thas. pont. 
1. Since the FFI is given only a Sinvte >sametemor ed 


zunction while the theoretically rous@creel francesa 
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looks at the complete function, the computer algo- 
rithm is designed to overcome this necessary short- 
coming by assuming the submitted finite sanpie 
repeats its pattern an infinite number of times. 
This can present probiems. If the right hand side of 
One pattern does not flow smoothly into the left hand 
side of the next repetition, as is shown in Figure 
3.4, then high frequency oscillations known as the 
Gibbs Phenomenon can be introduced into the spectrum. 
(See Figure 3.5) This is avery real danger where 
accustic pressure fields are concerned. For example: 
if the sample size doeS not include the entire 
Eaetern (1.¢., the pressure magnitude has not been 


completely attenuated prior to the endpoints of the 


Sample), this unwanted "jitter™ in the spectrum can 
result. TO | avotdmiinis Situation, not only was a 
pressure field symmetric about tne origin (tke 
source} used, but also a Hanning Window vas 


constructed to reduce both the amplitude and the 
phase components of the pressure field to zero at the 
endpoints, thus avoiding the Gibbs Phenomenon. 
Figure 3.6 shows the same pressure field depicted in 
Figure 3.4 combined with a Hanning Window. Its spec- 
trum is shown in Figure 3.7. Yhe Gibbs Phenomenon 1s 
absent. The form of the window is given below. The 
reader's attention is directed to Appendix I fora 
view of the sganner in which it is combined with the 


pressure field. 


The Hanning Window: 


WE 0.5| aaa coe (egn 3.1) 
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where "NPTS" ais the umber of data poirts, 


and [= Bienes. 


The Hanning Window does not affect the acoustic fres- 
sure field in any other way; the interference pattern 
in the range domain characteristic of the Lloyds 
Mirror for an Acoustic Doublet rewaains the sane 
except for its amplitude near the end points. AS 
Shall be shown later in this chapter, 1t is the 
pattern of the entire pressure field, not simply the 
amplitude of the pattern, oo which the ®.7T. performs 
its legerdemain. 

Because the vertical wavenumber beta is real and has 
physical meaning for values of the horizontal 
wavenumber less than k, only those values of the 
magnitude of the acoustic pressure spectrum -—sorre- 
Sponding to Y < k were retained. 

In the Wavenumber Technique, the relationship betyeen 
the horizontal range step Size A\c and ere hoki 
zontal wavenumber step size Ay is extremely 


important. Recalling equation 2.36, 





© 2Gi 
Ay &r N (eqn 3.2) 
for a given "N," the term on the right hand side, 
271 /N , is a constant. Conseguently, 
= 
rae es 17) (eqn 3-3) 


Simultaneously, in accordance with sampiing theory, 
in order to prevent aliasing in the pressure spec- 


trum, there iS an upper limit on the size of Ee 
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as 


(eqn 3.4) 
Normally, fA\r is chosen such that it is much less 
than half a wavelength; specifically, Ac < A/5 or 
/6, thus agood sample of the range domain is 
ensured. However, from eyuation 2.36, the smaller 
Ar becomes, the larger AY BeCole eal Guage 
coarser the sample sizes of the spectrum bécone. 
Llso, to ensure adequate samples in the spectrun, 
examination "62 equations 2225 and 2.36 and { Ref. 10} 


reveals that 


2 
Ay = sr (+) (egn 3.5) 


where Z is the deeper of the two depths, whether 
source or receiver. So there is a delicate balance 
to maintain between the sizes of A\r and Ay . en 
the Splbat Of Compromise, the investigator chose to 
use a range step size equal to a guarter of a wavel- 


ength where 


(egn 3.6) 


For most cases investigated, values for k,N, De 
and AY are Summarized in Table I. As can be seen, 
uSing these values and equation 3.5, for a "Z" value 
of 25 wavelengths, the maximum size AY can assune 
Gena stat ensure an adeguate sample size is 
Ay = 0.0004 /n. GWearly tie si zesior Ay shown 
mre tapes I are too coarse by a factor sof appEoxi- 
mately 20. Only hy Increasing ~N., the number of 


sample points, can this discrepancy be corrected. 


3.) 


B. 


However, aS mentioned earlier, computer resources 
were limited. The effect on the research results is 
not significant; nevertheless, it is recommended the 
reader bear this limitation in mind during the 
succeeding pagés. 

Looking again at Figures 3.1 threugi 3.27 ene ees 
the curve marked by the solid line, while the theo- 
retical Fourier Transform is the curve marked by the 
dotted line. Notice the large discrepancy between 
the theoretical values and tae actual vaiues in the 
region corresrfonding to small vertical wavenumbers. 
This happens to correspond to a region o£ the spec- 
trum characterized by more rapid fluctuations in the 
spectrua. However, the size of beta 1S consrane 
Conseguently, fewer samples of this region of the 
spectrum are available as compared with that region 
corresponding to high beta vaiues. Sincevthe. piet. 
ting routine utilized is interpolating between 
sampled points of the spectrum, the graph at this erd 
is less smooth and results in a marked difference 
between the curves. where the curves are determined 
by more points and are smoother, the discrepancy 
between theory and computerized reality is much less 


noticeable. 


SCURCE DEPTH DETERNINATION 


Theory predicts that, for a source Shallower than the 


receiver, when the magnitude of the pressure spectrum is 


plotted as a function of the vertical wavenumber, the resul- 


tant nodal spacing can be used to determine source depth by 


either consulting a chart (see Figure 3.8) or by performing 


a Simple calculation (see equation 2.31). To test this 


prediction, several scenarios were designed. Figwnes..3.9 


38 


Pieewgh 3. 1ieare “typical of the many cases run. In each, 
source depth is less than receiver depth. For each case, 
carerul measurement of the AB Spacing Goupled wath use of 
equation 2.31 yielded the (already known) source depth with 
Sesemethatea two per cent errorg® which *can be attributed to 


human error in the measurement of (Ny Spe 


TABLE I1 


Results of Source Depth Determination Runs 


Figure Nr. Zoey’ ZR (M) (Me als: ice aulecutclly 
3.9 ee, me 
3.10 CDS © ee 2 31.4" 

5 15703. kW ET) iS. 7 ™ 
3.12 47.124 21.998 o.14 " iit 
Sig 47.124 31.416 0.10 " See: 
ame allah 15.708 0.20 " 15.7" 


To test what happens when the receiver depth is shal- 
lower than the source depth, the cases depicted in Figures 
3.12 through 3.14 were run. The scenario is the same as for 
Figures 3.9 through 3.11; only the Source and receiver 
depths have been exchanged. AS can be seen, the AL 
Spacing corresponds with the receiver depth. This 1s again 


as predicted by theory. Table II summarizes tne findings of 


Pegiises 3.9 through 3.14. 


BS) 


The scenarios shown in Figures "S03 through Sei) a7 
start with a pressure field measured horizontally from the 
source; i.e., the exact range to the source is Known. What 
if the range to the source is not known, or oaly minimal 
information concerning the range is known? in, thas case, 
only a portion of the pressure field starting at scme 
Initial eance Po can be sampled. What 2s the e€£eecc men 
the spectrun? 

TO answer these questions, a range window was 
constructed and placed at varying initial ranges from the 
source. AS Figure 2.2 illustrates, as range from the source 
increases, the nodes of the Lloyds Mirror interference 
pattern in the range domain become more widely spaced. 
Given fewer nodes to sample, will the spectrum be the same? 
In accordance with theory, the spectrum is independent of 
range (see equation 2.39) and so should be the same wherever 
the window is placed. 

Figures 3.15 through 3.17 illustrate tiak the spectaam 
is not entirely independent of range. the figures utilize 
the Same scenario as was used above, differing only in the 
range from the source at which the sampling of the pressure 
field begins. AS can be seen, beginning at the right hand 
side which corresponds to high vertical wavenumber values, 
or the low horizontal wavenumber values (see eguation 2.8), 
a “washing out" of the spectrum cecetrs, increasing as the 
range window is moved further from the source. Again, the 
theoretical curve is marked by the dotted line. 

When the magnitude of the pressure spectrum is plotted 
aS a function of the vertical wavenumber, the useful portion 
of the resulting graph is the last two-thirds of the beta 
range. For K = 2 20 7ally this is the range 
QR lel yin < LB < 207s This phenomenon is explained more 
fully in Section III A 4. With this caveat in mind and 


based upon multiple source, receiver, Tange window 
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combinations run through the model, the following criterion 
waS eStablished: the nodal spacing of the spectral flots 
was no longer determinable once the initial range % was 
increased to approximately three times the source denth. In 
other words, minimum knowledge of the range from receiver to 
source must be availakle to ensure adequate sampling of the 
acoustic pressure field in the range domain. This, in turn, 
will produce a pressure spectrum of such a quality that 
source depth determination can be made. 

What causes this not entirely unexpected result? As the 
range window is moved further from the source, for 
Z. = Z , those samples of the pressure field corresponding 
to low gamma values are lost to the spectrum first. Low 
fauna values correspond to waveS arriving at high angles 
with respect to the horizontal. At greater receiver ranges, 
& more closely approximates r and the arrival angles beccocne 
@uoscer tc the horizontal (see Figure 2.1). 

As the window is moved further from the source and fewer 
nodes of the interference pattern are encountered, if a 
wider window were used, could this “washing out" erfect in 
the spectrum be minimized or even eliminated? ce 
aimenor’sS Opinion, it would be minimized; because of the 
reasoning in the preceding paragraph, it probably would not 
be eliminated. However, a wider window would require more 
sample points; because of equation 2.36, one Cannot merely 
make Ar larger. Time and the computer resources available 
Pemuiis investigator were limited, and the test cases run 
were already using the maximum availabie array sizes, tnus 
Pace igduna@ma  Guantitatiye illustration of this hypothesis. 
Succeeding research into the W.7T. might include ways to test 
times POLnt. 

In conclusion, the Wavenumber Technigue has the fpoten- 
tial for being a valuable operational tool in determining 


the depth of an acoustic source provided the receiver is at 
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a greater depth than the source. It is strongly reccnumended 
that further research into use of the W.T. in a realistic 


ocean environment be done. 


C. THE EFFECT OF SURFACE ROUGHNESS 


This and the succeeding section were purely gquaiitative 
portions of the reséarch. Therefore, in the figures the 
spectra were plotted as functions of the  hcrizoneat 
wavenumber for ease of discussion except where a particular 
point about the vertical wavenumber was being illustrated, 
and the theoretical curve was omitted. 

As deScribed in Chapter II, the effect of surface rough- 
ness is to suppress the contribution of the image to the 
dipole interference pattern. Moreover, by increasing M, the 
effect is to supvress the image's contribution to the pres- 
sure field in the range domain and reduce the pressure spec~ 
trum to the contribution of the direct path wave only (see 
equation 2.41); the resultant magnitude of the spectrum is 
inversely proportional to the verticai wavenumber (see egua- 
tlone2e42 As Figures 3.18 thzough 3921) 123ustrarcr eee 
is exactly what happens in the FFT. In Figure 3.2, am 
inverse beta curve has been manually superimposed onto the 
actual curve to reflect this point. 

Looking at the pressure” field in™=evesrange Gold meres 
the effect of surface roughness is inversely proportional to 
the wavelength of the sigzal, For longer wavelengths, the 
Lioyds Mirror effect is more tolerant of surface roughness 
than it is for shorter wavelengths. Does the same relaticn- 
ship hold for the fressure spectrunr? Comparison of the 
respective sea states in Figures 3.22 and 3.23 with Figures 
3.24 and 3.25 illustrate that it does. 

This investigator also looked at the effect of varying 


the range window in the presence of waves. This is similar 
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Dw idwmewas looked at an Section B of this chapter. Results 
should be similar and for the same reasons. For a given sea 
State, source and receiver geometry and wavelength, the 
range window was moved successiveiy further from the source 
(see Figures 3.26 through Sr 29) . AS can be seen, the 
“washing out" effect cf the spectral nodes begins to effect 


the lower horizontal wavenumber values first. 


tee 2h EFFECT “OE*NOISE 


Random noise was simulated in this research. Being 
independent of range and frequency, it was treated asa 
constant by both the Fourier Transform and the FFT. The 
investigator wished to compare qualitatively the effects of 
equal amcunts of noise on the pressure field in the range 
domain and in the wavenumber domain. TO do times ae CCcs- 
Sively larger values of the “noise factor" were used in 
computing n (see equation 2.43). FYalues ox the noise 
factor used in the research were limited to the maxinun 
amplitude cf the pressure field P(R). Since pol emrad? da 
t10n iS predicated, as range from the source increases, the 
maxima of the interference pattern decrease as the inverse 
EGuare of the range ({i1.e., oe Hence, even small noise 
factors can have a devastating effect on the pressure field 
aS the receiver iS stepped out in range. The destructive 
effect on the spectrum of sSucceSsively more intense amounts 
SuenmoOise at varying wavelengths is illustrated in Figtires 
S230 through 3.47 

SevenaimGcoucilUuSlons sconcerfing the pressure spectra’ Can 
be drawn from these results: 

1. For equal amounts of noise, longer wavelengths are 
affected less than shorter wavelengths. Compare 
Figures 3.37 and 3.32 with Figures 3.43 and 3.454, and 
Figures 3.34 and 3.35 with Fiyures 3.46 and 3.47. 
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This was rot a Surprising result since the same fprin- 
ciple holds for the range domain. 

2. The spectrum is affected less by noise thar is the 
range domain. Compare Figures 3630) tinoiciheS om 
Figures 3.36 through 3.47 for an illustratienser this 
point at two different wavelengths and various noise 
levels. This was a rather welcome surprise. 

3. If the maynitude of the pressure spectrum is plotted 
as a function of beta, the AL Spacsinagmcam still 
be detected even after tne pressure function in the 
range domain has been "swallowed up" oy the noise 
field. Compare Figures 3.33 tirewen. 3-35 eeneue eee 
3.42 through 3.44. This was not surprising in view 


of 2. above. 


E. SUMMATION 


All of the foregoing results are based specificaliyv cn 
the Lioyds Mirror for ar Acoustic Dewslet. Thess acea nee. 
idealized and artificiai environment. Any thougnt of imnme- 
diately applying conclusions reached in this paper to a real 
world situation is premature. Howevec, the results are 
Still significant if only for supporting the statement made 
by each preceding investigator that the Wavenumber Technigqne 
Should be examined very closely for a future operational 
role, especially in view of the current trends in scurce 
levels and ambient ncise levels. 

Certain fundamental conclusions regarding the author's 
research into the Wavenumber Technique can now be made 
within the confines of the above restriction: 

1. Source depth can be determined aguickly and easily 

xrom the acoustic pressure spectrum provided 


a) the receiver is deeper than the source, 
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b) some knowledge of the range from receiver to 
source is available,and 

©) the Magnitude of the pressure Spectzum is pletted 
as a function of the vertical wavenumber. 

Source depth determination in even a fhoisy environ- 

ment is possible. While the INtredwetrcn GE 

lncreasing amounts of narrow band noise adversely 

affects both the range domain andthe wavenumber 

domain, the pressure spectrum seems able to withstand 

the chaotic destruction longer than does the rressure 

field in the range domain. 

The ability to determine source depth is adversely 

affected by the height of surface gravity waves since 

surface roughness reduces the Lioyds Mirror as the 


sea state increases. 
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Range Window Set at 62.8 Meters 
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Figure 3.19 Pressure Spectrum, Sea State 2 
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Figure 3.25 Pressure Spectrun, Sea State 3, K = 2.0 
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Figure 3.27 Pressure Spectrun, 
SS 2, Range Window Set at 100 Meters 
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Figure 3.28 Pressure Spectrun, 
SS 2, Range Window Set at 200 Meters 
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Figure 3.34 Pressure Spectrum vs. Gamma 
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